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Einleitung. 

Tancht man einen Körper in eine strömende Flüssigkeit, so 
beobachtet man, daß vor dem Körper die Strömung glatt ist, hinter 
ihm dagegen die Flüssigkeit sich in starker Rotation befindet. 
Man kann sehen, daß diese Rotation nicht sofort eintritt, sobald 
sich der Körper relativ zur Flüssigkeit zu bewegen beginnt, 
sondern erst nach einer gewissen endlichen, wenn auch kurzen, 
Zeit. Bis zu dieser Zeit haben wir auch hinter dem Körper an- 
genäherte PotentisJströmung. 

Es ist schon vielfach versucht worden diese Erscheinung zu 
erklären und zu berechnen; doch gelang es erst in neuerer Zeit 
L. PrandtP) eine zufriedenstellende Theorie über die Entstehung 
der Wirbel zu finden. In kurzen Zügen sei diese angeführt *). An 
der Oberfläche des Körpers haftet die Flüssigkeit, hier ist also 
die Geschwindigkeit gleich Null. Mithin existiert eine gewisse 
Schicht — wir nennen sie Grenzschicht — , innerhalb der die 
Geschwindigkeit von Null auf den Wert der Potentialströmung an- 
wachsen muß. Ist die Zähigkeit der Flüssigkeit sehr klein, so ist 
auch die Grenzschicht sehr dünn. Wir haben in ihr also ein großes 
Geschwindigkeitsgefälle, was dazu führt, daß das Reibungsglied 
der hydrodynamischen Grundgleichung selbst noch bei kleiner 
Zähigkeit eine wesentliche Rolle spielt. 

Wenn die Grenzschicht in ein Gebiet höheren Druckes ein- 
tritt, so erfahren die Teilchen in jeder Entfernung von der Wand 
Verzögerungen; und da die Teilchen im Inneren der Grenzschicht 
von vornherein geringere Geschwindigkeit besitzen, so wird sich 
hier eher als draußen die Geschwindigkeit umkehren und damit 



1) Verh. d. intern, math. Kongr. 1904, p. 484 ff. 

2) Ausführlicher bei H. Blasius. Dissertation: Grenzschichten in Flüssig- 
keiten mit kleiner Reibung. 1908. Zeitschr. f. Math. u. Phys. p. 1 ff. 
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ZQ ßückstrSmnng and Strahlbildnng führen, wie die für verschie- 
dene Querschnitte gezeichneten Greschwindigkeitäprofile ^) in Fig. 1 
zeigen. Diejenige Stelle, die die Grrenze bildet zwischen den G-e- 



Fig. 1. 

schwindigkeitsknrven , die nnr positive, und denen, die auch ne- 
gative Teile haben, nennt man Ablösnngsstelle. In ihr muß das 
Oeschwindigkeitsprofil die Achse tangieren. Diese Bemerkung 
erlaubt nns, die Ablösongsstelle zq berechnen. 

Diese von L. Prandtl anfgesteilte Theorie warde von H. 
Blasins') besonders nach der mathematischen Seite hin weiter 
ansgebant. In der betreffenden Arbeit^) wird sie auf zweidimen- 
sionale Probleme angewendet, und es stellt sich dabei heraus, daß 
anch quantitativ mit Hülfe der oben genannten (rrenzschichten- 
theorie die Wirbelbildung hinter festen Körpern in stromenden 
Flüssigkeiten dargestellt werden kann. 

Die vorliegende Arbeit hat den Zweck , die Entwicklungen, 
die von H. Blasias für zylindrische Körper aufgestellt wurden, 
auf Botationskörper zn übertragen und gleichzeitig aaf mehr 
Glieder za berechnen. Wie später gezeigt werden wird, läßt 
sich ersteres nicht mehr auf analytischem Wege erreichen, wes- 
halb ich die Zuflucht zu numerischen Methoden nehmen mnß; zu 
deren Ansführnng wurde mir vom Institut für angewandte Mechanik 
zn Göttingen eine Rechenmaschine in freundlichster Weise zur 
Verfügung gestellt. Die Anregung zu dieser Arbeit verdanke 
ich Herrn Prof. Prandtl. Ebenso schulde ich Herrn Prof, Runge 
Dank, der mich bei den in dieser Arbeit verwendeten numerischen 
Methoden mit seinem Rat unterstützt hat. 



1) L, Prandtl 1. c. 

2) H. BUsiuft, DisMrUtioD 1908. (I.e.). 
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i. Formulierung des Problems und Aufstellung der 

Grundgleichungen. 

Wir nehmen die Konstante der inneren Reibung klein an. 
Wie bereits oben erwähnt, wird damit auch die Dicke der Grenz- 
schicht klein und zwar sei ihre Größenordnung mit s bezeichnet. 
Biß nahe an diese heran herrscht dann Potentialbewegung. 

1. Der Arbeit von H. Blasius gegenüber, soll in dieser 
das Problem der Rotationskörper, deren Achse in Richtung der 
ungestörten Strömung liegt, behandelt werden. Als Koordinaten 
für die Grenzschicht bieten sich als zweckmäßig die folgenden 
(Fig. 2) : 




Fig. 2. 

x = Bogenlänge längs einer Meridiankurve, 

y = normaler Abstand von der Oberfläche des Körpers; 
die Gestalt des Körpers selbst wird durch 

r = den senkrechten Abstand von der Rotationsachse, beschrieben. 

r wird in der Grenzschicht zunächst Funktion von x und y, wird 
aber schließlich nur in seiner Bedeutung als Querschnitt- 
radius des Rotationskörpers (y = 0) verwandt. 

2. Auf diese Koordinaten müssen wir die für rechtwinklige 
Carthesische Koordinaten aufgestellten hydrodynamischen Grund- 
gleichungen, die in der Vectorschreib weise lauten: 

öö . , _, 1 ^ Je ^ 
-^r + (l'r)ö = rp-^ — Jt) 

Ot ^ Q Q 

div = 
umformen. Die Umformung der Beschleunigungsglieder (Dr')ö 

liefert neben Gliedern wie u -^ , v _- u. s. w. Zentrifugalkräfte 

des Typus — '— und — ; diejenigen von ihnen, die in der Gleichung 
für die a;-Komponente auftreten, dürfen, da sie von der Größen- 
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Ordnung s sind, gegen u — + v — r\jl fortgelassen werden können. 

Die Resultierende der Reibungsschubkräfte ist nach neben- 
stehender Figur 3 
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Aus der dynamischen Grieichung für die a;-Richtüng bleibt also 
übrig: 

du , du du 1 dp k d*w 

dt dx dy Q dx Q dy* 

Die y-Gleichung liefert bei den Größenordnungsabschätzungen 

-^— r\^ — ro 1. Hieraus schließen wir weiter , daß -:r- i^ der 
dy r ' dx 

Grenzschicht von -^ in der Potentialströmung nur um eine Größe 

von der Ordnung $ verschieden ist, die wir in der a?-Gleichung, 
wie wir ja nur Glieder r^ 1 beibehielten, vernachlässigen. Somit 

dj) 
können wir ^ in der Gleichung aus der Potentialströmung be- 
rechnen, also für die Grenzschicht als eingeprägt betrachten. Die 
y-Gleichung ist damit für unsere Zwecke vollständig ausgenutzt, 
sie würde anderenfalls nur zur subtileren Berechnung von p in 
der Grenzschicht dienen. 

3. Die Kontinuitätsgleichung wollen wir schreiben 

du dv ^^ __ (\ 

wobei wir für einen Augenblick 5, i?, E als die rechtwinkligen 
Koordinaten betrachten. Ganz allgemein nehmen wir dann die 
Transformation: 



i 



) 
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vor. Durch ganz elementare Beclinnng findet man die neue Kon- 
tintiitätsgleicliang : 



dx \ dxyz 
nnd zwar bedeutet: 



dxyz 



u. 



dx 



dx dy dz 

drj drj drj 
dx dy de 

dx dy dß « 
dy r de 

dt' ^' dt 



Die gefundene Gleichung ist vollständig allgemein gültig. 

Auf unsere Koordinaten der Rotationskörper angewendet gibt 
dies folgendes: 

q> ist Azimut (Fig. 4). 
i 




Fig. 4. 



Wir haben also: 



V = r(x,y)Bm<p 
t = ^(^3 y)cosg) 

Berücksichtigt man, daß x, y und q> senkrecht auf einander stehen, 
so ist der Wert der Funktionaldeterminante: 
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dxytp 



= r. 



Da femer ans Syznmetriegründen keine Strömnng normal znr xy- 
Ebene stattfinden kann, so muß unsere Kontinuitätsgleichung 
lauten : 

In unserem speziellen Falle hätte man dies Resultat auch direkt 
durch geometrische Überlegung finden können, indem 

-^ (2^r dy u) + — (2nr dx v) ^ 

sein muß. 

Ist die Dicke der Grrenzschicht von der Grrößenordnung s, dann 
ist 



nu 



1, -^ '^ 1> Vf, r\J e. 



dx dy 

Auf die Kontinuitätsgleichung angewendet: 

du^ dr dvo , dr 

2 2 2 2 2 2 2 

11 11 11 8 

Da — als Cosinus der Neigung der Meridiankurve von der 

Größenordnung! ist, darf das letzte Glied vernachlässigt werden; 
r ist mithin nur als Funktion von x zu betrachten. 

4. Bezeichnen wir nun auch in der Kontinuitätsgleichung die 
Geschwindigkeit mit w, v, so haben wir als Differentialgleichungen 
der Strömung um Rotationskörper die folgenden: 

fdu, du du\ (du dil\ , ö"tt 



( du , du du\ (du dü\ , 



dy' 



dru drv ^ 



dx dy 

mit den Grenzbedingungen : 

für y= 0:w = 0, f; = 
für y = oo: u = ü. 

ü ist die Geschwindigkeit der Fotentialbewegung. 

Kennen wir die Integrale dieser Gleichungen, so beherrschen 
wir damit die gesamte Strömung um alle Rotationskörper, sobald 
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uns die Potentialströmxmg um diesen Körper bekannt ist. Zur 
Losung mnß gegeben sein, wie Körper nnd Flüssigkeit sich gegen- 
einander bewegen ; also entweder die Geschwindigkeit, mit der der 
Körper durch das Wasser gezogen wird, oder die Strömungs- 
geschwindigkeit des Wassers, in die der Körper getaucht wird. 
Beides ist mathematisch nur wenig verschieden. 



IL Entstehung der Grenzschicht und Ablösungsstelle beim 
plötzlichen Beginn der Bewegung aus der Ruhe. 

1. Wir wollen das einfachste Bewegungsproblem wählen. 
Der Körper ruhe bis zur Zeit ^ = in der Flüssigkeit. Zur Zeit 
^ = erhalte er plötzlich eine bestimmte Geschwindigkeit ü und 
werde dann ohne Beschleunigung weiter bewegt. Die Schwierigkeit 
der Ausführung des plötzlichen Stoßes kann man dadurch umgehen, 
daß man den Körper nicht plötzlich, wohl aber hinreichend schnell 
auf die Geschwindigkeit ü bringt, d. h. so schnell, daß die Zeit r, 
die dazu erforderlich ist (Fig. 5), klein gegen die Abi ösungs zeit 
(s. w. u.) ist. In diesem Falle sind nur geringe Abweichungen 
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Fig. 5. 

der Erscheinungen von den theoretischen Betrachtungen zu er- 
warten. 

Unter der Einwirkung des Stoßes stellt sich im ersten Moment 
Potentialbewegung ein mit Ausnahme einer unendlich dünnen Schicht 
an der Körperoberfläche. Von hier aus setzt unter der Wirkung 
der Reibang die Bildung der Grenzschicht ein. 

2. Die Gleichungen, die zu benutzen sind, lauten: 

du du du du ö'm 

dt dx dy dx dy^ 

dru drv ^ 

dx dy ' 
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h 
wobei X = — gesetzt ist. Die zweite Grleichtmg integrieren wir 

dnrch 

— iii* — 1 a» 

r dy r dx 

Die Grleiehnngen stimmen bis auf den Faktor — mit den Glei- 
chungen überein, die Blasius (1. c.) in dem entsprechenden Problem 
für den Zylinder fand. Wir können daher den dort gefundenen 
Ansatz entsprechend übernehmen, ohne die Überlegungen, die dazu 
führen, zu wiederholen. 

3. Wir machen den Ansatz 

y 



V = 



2\/xt 



Durch Einfuhren dieser Werte in die dynamische Grleichung findet 
man die Reduktionsformel:* 

»• ^ \^n\ äxdrj r dx drj ) dx dif J' 

für (1 = 1 tritt zur rechten Seite — 4wr— additiv hinzu. Geben 

wir (Jt die Werte 0, 1. 2 ... , so erhalten wir die Differential- 
gleichungen für X' 

da; r dij \öa;di? r da; öij / r öx dif 

Beschränken wir nns znnächst auf diese beiden. Jedes %, berech- 
net sich aus den vorhergehenden %, bis jj^,, sodaß die Gleichungen 
successive integrierbar sind. Da x, Funktion von x und y ist, 
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Tmsere Differentialgleichangen aber nur die Abhängigkeit des x, 
von 11 bestimmt, so innß sich %^ anfbaaen lassen ans einer Sninme 
von Produkten ans Funktionen von i; nnd Funktionen von x. 
Unter Rücksichtnahme anf die Grenzbedingong setzen wir: 

wo S, allein von ri abhängt und der Gleichung genügt: 

C + 2^;r. = 0. 

Die oberen Striche bedeuten Differentiation nach 17. 

Dies führen wir in die Differentialgleichung von %^ ein. Man 
findet, daß dann 

zu setzen ist, wo ^, den Gleichungen zu genügen hat: 

■c+2i?c-4s;. = 4(c-g.r.-i) =4A. 
5r;+2i?e;'»-4g;» == -45,5; = 4a». 

Aus genannten G-leichungen sind dann ^0, l^^y tih ™^^ ihren Ablei- 
tungen zu berechnen, und man findet dann 



t* = ^^i[ + t 



(«'S K.+»'l •£«.)+['■)■ 



Also haben wir nun die Aufgabe, folgende Grleichungen zu lösen: 

C + 2ijC =0 

c+2n;'.-4e;„ =4(c-&.g;-i) 
5;';+2i?G-4r.» = -4s.r; 

mit den Grenzbedingungen: 

für 1? = 0: g, = 0; g: = 0; t,. = 0; i[^ = 0; g., = 0; f ;, = 0, 

füriy = oo: ?; = i; e;. = 0; g;, = o. 

4. Von H. Blas ins (1. c.) sind die Gleichungen in geschlos- 
senen Ausdrücken für g^ und g,, gelöst worden, indem es möglich 
war, die homogenen Gleichungen zu lösen. Da wir dies Resultat 
später noch brauchen werden, möge es hier angegeben werden. 
Die Gleichungen für tn sind vom Typus: 

JF:; + 27^Fi-4nF, = 0. 

Die beiden Integrale sind dann durch folgende Ausdrücke be- 
stimmt: 
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Wie schon diese Ausdrücke zeigen, werden die g, je größer der 
Index wird, desto komplizierter und daher für weitere Becbnmigen 
unpraktischer. Dazu kommt noch die Schwierigkeit, daß bei » > 2 
die eingeprägten Kräfte in den Gleichungen für g^ Glieder ent- 
halten, die mehr als dritten Grades in e^^ sind. Diese lassen 
sich nicht mehr zurückführen auf die Funktion e ^ und 

CO 

multipliziert mit. Potenzen von t], was zur weiteren Rechnung er- 
forderlich sein würde. Es bot sich daher ein numerisches Ver- 
fahren als das nächstliegende dar, wodurch die zweite Schwierig- 
keit sicher gehoben wird, und beliebig viel Glieder in der Ent- 
wicklung von t sich berechnen lassen. 

5. Die Methode, die ich zur Berechnung verwendet habe, ar- 
beitet mit den verschiedensten Hülfsmitteln. Die Gleichung 

C + 2nC = 

kann man exakt integrieren: 

woraus durch einmalige Integration folgt: 

^0 = ß + ^f e "^ dn. 

oo 

Die Konstanten bestimmen sich aus den Grenzbedingungen 

für 71 = 0: to = 
für iy = OO: g*^ = 1 



Also: 






\ 
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sodaß : 

Die Funktion berechne ich punktweise und leite daraus die In- 
tegrale Co ^°^ 5o durch numerische Integration ab. In diesem 
FaUe, wo ^^ noch durch eine verhältnismäßig einfache Funktion 
von ri dargestellt ist, würde es nicht weniger einfach sein, das- 
selbe aus der Formel 



f. = i+>/V'V, 



zu berechnen. Jedoch schon für ^ wird sich die numerische Rech- 
nung einfacher gestalten. 

6. Die numerische Integration ist nach der Simpsonschen 
Regel ausgeführt. Wir setzen 






Dann ist angenähert 

Dieser Nährungswert weicht vom wahren Wert des Integrales um 
Unterschiede ab, die mit l zugleich unendlich klein werden von 
der vierten Ordnung. Bei dem zu behandelnden Problem sind die 
Kurven berechnet von 0,1 zu 0,1 Einheiten der 17 Koordinate. 
Es ist also: 

l = 0,05 

l' = 6,26.10-. 

Da die Kurven nicht stark oscillierend sind , so wird uns die 
Simpsonsche Regel die Werte der Integrale mit genügender Ge- 
nauigkeit geben. Die Integrationskonstanten sind immer durch 
die Grenzbedingungen bestimmbar. So sind nun die drei Kurven 
S'o? So- to numerisch berechnet. Die Werte sind in der folgenden 
Tabelle ^) in Intervallen von 0,1 zu 0,1 Einheiten der iy Koordinate 
angegeben (Fig. 6). 

1) Die Werte der Integrale fe^'^dx sind entnommen : A. Markoff, Tables 

OD 

des valeurs de Fint^gral T e~'*<Ä. 

m 
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r» 




V 







0,1 



0,2 0,3 



0,4 



0,5 



0,6 



t: 


1,128 


1,112 


1,084 


1,031 


0,962 


0,879 


g, 





0,112 


0,222 


0,326 


0,427 


0,519 


to 





0,006 


0,022 


0,050 


0,087 


0,136 



0,787 
0,603 
0,191 



0,7 



0,8 



0,9 



1,0 



1,1 



1,2 



1,3 





0,691 


0,596 


0,502 


0,415 


0,336 


0,267 


s; 


0,676 


0,741 


0,796 


0,842 


0,879 


0,910 


So 


0,256 


0,326 


0,403 


0,484 


0,571 


0,660 



0,208 
0,934 
0,752 



V 



1,4 



1,5 1,6 



1,7 



1,8 



1,9 



2,0 



r; 


0,159 


0,119 


0,087 


0,063 


0,044 


0,031 


r. 


0,952 


0,966 


0,976 


0,984 


0,989 


0,993 


So 


0,846 


0,942 


1,039 


1,137 


1,236 


1,335 



0,021 
0,995 
1,434 



7. Mit Hülfe dieser Zahlenwerte können wir die eingeprägten 
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Kräfte der Differentialgleichrmgen fiir tia ^^^ tib bereclmen. 
war: 

4f„ = 4(C-SX-1) 

4U = -45.C 

Wir erhalten die folgenden Werte (E^. 7 u. 8). 



Es 




Fig. 7. 



«J 





0,1 


0,2 


0,3 


0,4 


0,5 


0,6 




4,000 



3,975 
— 0,025 


3,850 
— 0,097 


3,776 
— 0,205: 


3,607 
— 0,337 


3,396 
0,474 


— 3,149 
0,601 


12 


0,7 


0,8 


0,9 


1,0 


1,1 


1,2 


1,3 


4/-„ 
4A» 


2,874 
0,705 


-2,584 
— 0,775 


— 2,273 

— 0,808 


1,969 
0,804 


— 1,673 
0,768 


1,395 
-0,706 


1,141 
0,626 


»J 


1,4 


1,5 


1,6 


1,7 


1,8 


1,9 


2,0 


4/ 


- 0,914 

— 0,638 


0,717 
0,448 


0,552 
0,363 


— 0,416 
0,285 


— 0,307 
0,218 


0,222 
0,163 


0,157 
0,119 



8. Wir können jetzt znr Integration der Differentialglei- 
chungen schreiten: 



1 
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mit den Grenzbedmgangen : 

füriy = 0: g,. = 0; t[a = 0; g,, = 0; i[, 

für 1? = oo: g;^ = 0; t[. 




0. 




Fig. 8. 

Die Integration kann mittels des Kuttaschen ^) Verfahrens aus- 
geführt werden, dessen Schema kurz angeführt werden soll. 
Der Grleichungstypus ist der folgende: 

Diese Differentialgleichung zerlege ich in zwei simultane Differen- 
tialgleichungen. Ich setze 



dK 
dfi 



y 



führe also y als neue Variable ein. Dann wird 

dV "" dri ' 
sodaß wir folgendes ßleichungssystem haben: 



dF, 
dri 



= y- 



1) Zeitachr. f. Math. a. Phys. Bd. 46. p. 436 ff. 



. 
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Wir schreiben zur Abkürzung: 

y = Y,iv,Y„,y). 
Ich bilde sodann nach einander die folgenden Ausdrücke: 

h = r,(iy, F„, y)h 

K = Yxiv + h, i^n+Ä,, y+^A 

Diese Ausdrücke kombiniert man zur Bildung neuer in der Weise, 
daß man setzt: 

Pt = kiK+K) Q. = HK+h) 

* = Pf+h(jPi-Pt) ^ = g, + i ((?! - ff.)- 
Bei diesem Resultat sind k und l die Korrektionen, welche 

an F^ resp. , " angebracht werden müssen, wenn man iy um die 

Strecke h wachsen läßt. Die G-enauigkeit dieses Verfahrens ist 

die gleiche wie bei der Simponschen Kegel. Die geometrische 

Darstellung genannten Verfahrens ist durch Fig. 9 gegeben. Zur 

Ausführung dieses Verfahrens ist nur erforderlich, daß man in 

dem Punkte, in welchem man anschließen will, die Werte von F^ 

dF 
und , ** kennt. In unserem Problem beginnen wir die Kurven 

vom Nullpunkt, und in diesem sind uns in der Tat auch die beiden 

Werte: F„(0) und (-^] bekannt. 

Es ist offenbar, daß dies Verfahren die Differentialgleichung 
um so leichter integriert, je größer der Krümmungsradius der 
Kurve an der zu berechnenden Stelle ist. Von der Kurve g( 

2 
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wissen wir nun, daß sie für i^ = und 1^=00 Null ist, also im 
einfaclisten Falle ein Maximum hat. Dann muß die £i Kurve, die 
für iy = ebenfalls gleich Null ist, parabolisch aus dem Nullpunkt 
herausgehen und sich dann einer endlichen Asymptote nähern; jeden- 
falls wird, wie das Kurvenbild später auch zeigen wird, die 
g Kurve flacher sein, als die g' Kurve. Wir werden daher ver- 
suchen, ob wir die Differentialgleichung nicht noch einmal inte- 
grieren können. In der Tat ist dies leicht möglich: 





Diese Form der Gleichung hat den weiteren Vorzug, daß sie nur 
die Ableitxmgen von g enthält, die man zu späteren Rechnungen 
braucht. 

Die Integrationskonstanten ergeben sich sofort, wenn man 
1^ = setzt: 

Cin = 0) = c,„- CC»? = 0) = c,„ 

5ia(0) nnd t[b(9) wären noch zu bestimmen. 




9. Die Berechnung dieser Größen gelingt mit Hülfe des 
Greenschen Verfahrens. Wir gehen aus von der ursprünglichen 
Gleichung : 

Die eingeprägte Kraft 4/,^ sehen wir als Superposition vieler auf 
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Strecken d^ wirkender Kräfte (Stöße) an ; dies können wir wegen 
der Linearität der Gleichnngen. Der Einfloß, den jeder einzelne 
Stoß z. B. anf t[a aasübt, ist proportional der Intensität des Stoßes 
und der Strecke di^, femer proportional einer Funktion d", die 
einmal abhängt von der Stelle ^, wo der Stoß stattfindet und von 
der Stelle rj, für die die Wirkung des Stoßes zu berechnen ist. 
Es wird dann: 



sein, d. h. gleich dem Integral über die Einflüsse aller Stöße. Man 
sieht leicht ein, daß 0* in 17 Lösung der homogenen Gleichung sein 
muß, da ein Stoß 4/',o(^)^^ nur auf eine plötzliche Aenderung der 
Integrationskonstanten hinwirken kann. Genau, wie ich d' so 
wählen kann, daß es mir Sla(^) bestimmt, so läßt es sich augen- 
scheinlich auch so wählen, daß g'/^ bestimmt wird. Wir suchen 
aber gerade g',a(i? = 0). Wir setzen also: 



Durch partielle Integration erhält man : 



ö,; 



Damit rechte und linke Seite der Gleichung identisch werden, muß 
d" der adjangierten Differentialgleichung 

genügen mit den Grenzbedingungen: 

für 1^ = : -ö- = — 1 

für t; = oo: #' = 0, ^ endlich. 

Nun sind die Lösungen der adjungierten Gleichung bestimmt durch 
die Lösungen der homogenen ursprünglichen Gleichung nach den 
Formeln : 

F F 

^ -^ FF' - F F' FF' —F F' ^ 

F^y^ und JP,a sind die oben angeführten Lösungen der homogenen 
DMFerentialgleichung (s. S. 12 0.). Die Determinante, die im Nenner 
steht — bezeichnen wir sie mit D — , genügt der Differential- 
gleichung : 

2* 
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D' + 2r]D = 0. 



Also 



D = e 



-1* 



Für d ergibt sich folgende Lösung: 

,1 r^ _«« 



» = tt 



{l + 2t}^e'^ f e-'' dri + Ti 



oo 



Das zweite Integral kommt nicht in Betracht, da d" im Unend- 
lichen endlich bleiben maß. Die Integrationskonstante a wird 
durch die Grenzbedingnng ^(0) = — 1 zu 

2 



a = 



Vä 



bestimmt, so daß d'y das znm Unterschied von später za berech- 
nenden mit «d-io bezeichnet werden soll, jetzt lautet 



YÄ L cx) J 

Mit Hülfe dieser Punktion können wir nun sofort ^[^ (0) berechnen. 
Es ist 





In ganz analoger Weise ergibt sich für gij(0): 



Führen wir diese Rechnungen aus, so erhalten wir folgende Ta- 
belle (Fig. 6, 7 u. 8) : 



V 





0,1 


0,2 


0,3 


0,4 


0,6 


0,6 




1 
4 



0,802 
3,186 
0,020 


0,648 
2,527 
0,063 


— 0,528 
1,995 
0,108 


0,434 
1,665 
0,146 


0,369 
1,220 
0,170 


— 0,300 
0,943 
0,180 


V 


0,7 


0,8 


0,9 


1,0 


1,1 


1,2 


1,3 




0,251 
0,723 
0,177 


— 0,212 
0,549 
0,166 


0,181 
0,410 
0,146 


— 0,154 
0,304 
0,124 


— 0,133 
0,222 
0,102 


— 0,116 
0,160 
0,087 


— 0,100 
0,114 
0,062 
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V 


1,4 


1,5 


1,6 


1,7 


1,8 


1,9 


2,0 


4/;„*.« 

4/",»*,o 


0,087 
0,079 
0,047 


0,07»i 
0,055 
0,034 


-0,067 
0,037 
0,024 


0,059 
' 0,025 

0,017 

1 


— 0,052 
0,016 
0,011 


-0,047 
0,010 
0,008 


— 0,042 
0,007 
0,005 



Durch Integration erhalten wir sodann: 

C{0) = 1,614 
C (0) = 0,169. 

Mittels des Greenschen Verfahrens können wir weitere Q-renz- 
bedingongen aasrechnen und zwar die Asymptoten, welchen sich 
die Kurven g,a und g,j im Unendlichen anschließen. Wir setzen 
dem obigen entsprechend: 







..oo 





Da die Grenzenglieder ^i^ noch nicht enthalten, integrieren wir 
noch einmal partiell: 



00 

1 

'o 






09 








Wir liaben also für &^^ die adjimgierte Differentialgleicliang : 



*;';-2i?»L-8»L =0 



loo 



mit den Grenzbedingungen : 

für 1? = 0: ^,. = 

für ./ = cxd: ^L-2)?^;„-6^,„ = 1; ^lo. endüch. 
Um diese Gleichung zu lösen, gehen wir von der Gleichung 

aas, deren Lösung 

*.« = - i + a 1? + (1 + 2,;') e'^'fe- "'d J 

oo J 
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ist, and es ergibt sich aus der G-irenzbedingang 

*„(0) = 0: 



sodaß 



BSJtc 



^i« = -i- 



3\Jn; 







Dies muß aber auch die Losung der adjungierten Differentialglei- 
chung 

sein, da man diese erhält, wenn man 

einmal nach rj differentiiert. -Ö-jo und -9",^ unterscheiden sich offenbar 
nur durch einige Eonstanten. Es ist 

sodaß man sehr leicht d'^^ aus ^,o berechnen kann. Wir haben 
damit 

/»CO 
5laW = 00)==./ ^loD '¥lJV 



und ganz entsprechend: 



tib iv = 00) =J -Ö-,« . 4ir,j dt], 



Diese Werte sind zur Lösung unserer Differentialgleichung von g^^ 
und g,, nicht erforderlich, wohl aber geben sie eine willkommene 
Eontrolle der gefundenen Eurven und haben als solche einen 
wohlberechtigten Wert. Die numerische Berechnung gibt folgende 
Werte: (Fig. 7 u. 8). 



V 





0,1 


0,2 


0,3 


0,4 


0,5 


0,6 







0. 


- 0,033 
0,131 
0,001 


0,0B9 
0,229 
0,006 


0,079 
0,297 
0,016 


0,094 
0,340 
0,032 


0,107 
0,363 
0,051 


— 0,117 
0,368 
0,070 



^' 
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'n 



0,7 



0,8 



0,9 



1,0 



1,1 



1,2 



».. 


0,125 


- 0,131 


0,137 


0,141 


0,145 


— 0,147 


4/;a*.« 


0,359 


0,339 


0,310 


0,277 


0,242 


0,206 


4/;»*,« 


0,088 


0,102 


0,110 


0,113 


0,111 


0,104 



1,3 

0,150 
0,179 
0,09 8 



*; 



1,4 



1,5 



1,6 



1,7 



1,8 



1,9 



2,0 



^.. 


0,152 


— 0,154 


— 0,156 


0,157 


0,158 


~ 0,159 


4/",„*.. 


0,139 


0,110 


0,086 


0,065 


0,048 


0,036 


4/'.j*.. 


0,082 


0,069 


0,056 


0,045 


0,035 


0,026 



0,160 
0,025 
0,019 



Durch Integration erhalten wir dann: 

eia(^ = OO) =r^^,^.4:f,aän = 0,413 



/oo 
».«•4/-.»rfi2 = 0,123. 



Wir kennen jetzt von unseren sechs Kurven ^[,, §;,, g,., g'^j, §;», g,j 
teUs dnrch gegebene Grenzbedingongen , teUs durch mit dem 
G-reenschen Verfahren berechnete, alle B.andwerte im Punkte 
ij = und im Unendlichen: 



für ij = 0: C = 1,614 

g.. = 
für ij = c»: g;', = 

r.. = 



i\, = 0,169; 

g;* = o 

g.» = 0,123. 



5.. = 0,413 

10. Jetzt können wir zur Integration der Differentialglei- 
chungen sckreiten: 

eL + 2i?L-6?,, =/ 4/;,diy + 1,614 = 4f,, 







.oo 



g^i+a'jr.j-eg,» = / 4/;,d, + 0,169 = 4f.». 



Die Werte von f,. und f,^ sind leicht zu finden und sollen daher 
nicht 'besonders angegeben werden. Wir wenden nun auf die 
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Differentialgleichnngen das oben angegebene Knttasche Verfahren 
an, wobei wir wieder 

Ä = 0,1 

annehmen. Wir erhalten so die Werte von gj und g'j. 

11. Es gilt nun aus i\a und g',j durch Differentiation die 
zweiten Ableitungen zu finden. Die numerische Differentiation 
wurde mit Hülfe der Differenzenrechnung ausgeführt. Auch diese 
Methode^) sei kurz angeführt. 

Die Kurve g(iy) sei als Punktion von tj in den Punkten 0, 0,1, 
0,2 u. s. f. bekannt. Will man uun für g eine angenäherte Kurve 
finden, die differentiiert werden kann, so schreibe man für den Punkt, 
für den man die Ableitung haben will, die folgende Entwicklung an : 

g = ao + a,A + |^Ä(Ä-0,l) + |^(A + 0,l)Ä(Ä-0,l) 

+ |^(A + 0,l)A(Ä-0,l)(Ä-0,2) 

+ |5.(Ä+o,2)(A + 0,l)A(A-0,l)(A-0,2) 
_i .... 

h ist dabei Abszisse, wenn man das Koordinatensystem in den 
Punkt legt, nm den man entwickelt. Die Koeffizienten a, sind 
vorläafig noch anbekannt, sie lassen sich aber leicht darch Diffe- 
renzenbüdnng bestimmen. Bezeichnen wir: 

^e = gM-g(A-o,i) 

dt = g(A+0,l)-5(ij„). 

Der Index bezeichnet dabei den Wert von ij, in dem der neue 
Anfangsponkt liegt. Man erhält dann: 

Jt = a.+|^0,2(A-0,l)+-^0,3A(A-0,I) + |i-0,4Ä(A-0,l)(A-0,2) 

+ |?0,5(Ä + 0,l)Ä(A-0,l)(Ä-0,2)+... 
SJt= 0,1. a, +o,l|i.O,2A +0,l|iO,3A(A-0,l) 



+ 0,l|?-0,4(A+0,l)Ä(A-0,l) + - 



^*^e = +(0,l)*a, +(0,l)'|i-0,2(Ä-0,l) 



1) C. Runge, Nature, Bd. 60 p. 865. 
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a 



+ (0,1)*-J|0,3Ä(Ä-0,1) + "- 



a. 



dJdJt= +(0,l)*a, + (0,l)« 210,2*+ ••• 

JSJdJi t= + (0, 1)« a. + • • . 

Daraus ergeben sich sofort die Koeffizienten zu: 

«0 = 5(Ä = 0) = m, 

a. = Jl{h = 0,1) = m, 

0,1 o, = dJlQi = 0) = «I, 

(0,1)' a, = JSJl{h = 0,1) = m. 

(0,1)' a^ = 8JSJt{h = 0) = •»»« 

(0,1)* a, = JSJdJiQi = 0,1) = m,. 

Diese Koeffizienten sind dann leicht ans folgendem Differenzen - 
Schema zn entnehmen : 

g(A = -0,3) g(A = -0,2) g(A = -0,l) t(h = 0) g(A = 0,1) g(Ä = 0,2) e(A = 0,3) 

V. 
^g(-0,2) ^5(-0,l) Jt(0) Jt{0,l) ^g(0,2) ^g(0,3) 

/ 
dJti-0,B) dJ^i-0,2) SJt(-0,l) dJtjO) <J^g(0,l) <y^S(0,2) 

\ 
^d^g(-0,2) JSJti- 0,1) JäJtiO) ^d^g(0,l ) -^d^f(0,2) 

/ 

«y^<j-^g(-o,2) d^d<^e(-o,i) ijijtjo) dJdJt{o,i) 

\ 

• • • • 

Die durch Unterstreichen hervorgehobenen Glieder sind unsere 
Koeffizienten. Nun ist es auch leicht, von g die Ableitung zu 
bilden. Dazu ist nur nötig, f nach h zu differentiieren. Im Punkte 
A = haben wir also: 

+ |^(0,2)(0,l)(-0,l)(-0,2) + ... 

0,1 (0,1)' . (0,1)'. 0,2 

^ 21 ■ 31 »^ 4! * 



ö! 
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»w, w- 2m. 4m. 
' 2! 3! ^ 4! ^ 51 ^^ 

Die Koeffizienten m sind dann genau die hervorgehobenen Diffe- 
renzen obigen Schemas. Unsere Rechnung verlangt nun nicht die 
Ableitung von ^ sondern von ^', die aber genau nach demselben 
Schema ausgeführt wird. 

Die letztgenannte Gleichung kann man noch etwas symme- 
trischer gestalten, indem man nicht als erste Differenz die Ope- 
ration J sondern d nimmt. Die Koeffizienten nehmen dann die 
Form an: 

< = g(A = 0) 

a[ = dtih = -0,1) 
0,1a; = Jd^{h = 0) 
(0,1)' oi = SJätih = -t),l) 
(0,1)' < = JSJ8i{k = 0) 
(0,1)* o; = SJÖJd^ih = -0,1) = m[ 

Man sielit, in dem Differenzenschema ist der Koeffizientenzng dem 
yorber gefandenen symmetrisch. Es ist dann noch: 





fw; 




m^ 








m[ 






^ 


^m 


— 


w?,' 




w. 






— ^ 


m^ 






' — 


JäJt{0) 




m[ 


— ' 


m^ 








♦Wj 






— 


JSJäJtiO). 



( 



*L\ - m; m;, 2m[ 4m; 

dAA = o ^ ' 2! 3! "^ 4! "^ 5 



Um die beiden Ausdrücke für f-rf ) zu einem zu verschmelzen, 
addieren wir diese Gleichung zu der obigen und dividieren durch 2 : 



Ua/» = o~ ' 2! 31'^ 4! ^ 



4^. . 

■* I ■ • • « 



51 

wobei 

_ m, + ml 

ist. Durch diese Operation bekommen die Koeffizienten im obigen 
Schema die folgende symmetrische Lage: 

g(-o,i) m g(+o,i) 

m[ Ml f»! 

wi Jtf, m, 
t«^ = Jl/^ 

tn; M, m,. 

Mit Hülfe genannter Methode kann man den Differentialquotient 
einer Kurve* mit der überhaupt erreichbaren Genauigkeit finden. 
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Die ersten Pankte von g", d. h. für 7/ = 0,1, 0,2 . . . die sich 
niclit mit dieser Methode finden lassen, falls man für negative 7; 
die Werte von J;' nicht kennt, kann man durch Reihenentwicklung 
um 1^ = leicht finden, da man ofi^enbar alle höheren Ableitungen 
von i im Punkte iy = berechnen kann. 

Man kann nun in Vollständigkeit die g^ -Kurven aufstellen 

(Fig. 10) : 



n 


&,. 


g;. 


SU 


tit 


bii 













1,614 








0,169 


0,1 


0,007 


0,142 


1,230 


0,001 


0,017 


0,170 


0,2 


0,027 


0,246 


0,872 


0,003 


0,034 


0,170 


0,3 


0,055 


0,318 


0,573 


0,008 


0,051 


0,164 


0,4 


0,089 


0,362 


0,310 


0,013 


0,066 


0,149 


0,5 


0,126 


0,382 


0,092 


0,021 


0,080 


0,126 


0,6 


0,164 


0.382 


0,081 


0,029 


0,091 


0,094 


0,7 


0,201 


0,367 


0,213 


0,039 


0,099 


0,057 


0,8 


0,237 


0,340 


0,305 


0,049 


0,103 


0,017 


0,9 


0,269 


0,307 


— 0,360 


0,059 


0,102 


0,021 


1,0 


0,298 


0,269 


0,384 


0,069 


0,099 


- 0,055 


1,1 


0,323 


0,231 


— 0,396 


0,079 


0,092 


0,081 


1,2 


0,344 


0,191 


0,367 


0,087 


0,083 


-0,100 


1,3 


0,362 


0.158 

* 


0,335 


0,095 


0,072 


0,109 


1,4 


0,376 


0,126 


0,298 


0,102 


0,061 


0,111 


1,5 


0,387 


0,099 


0,253 


0,107 


0,050 


0,107 


1,6 


0,396 


0,075 


0,210 


0,112 


0,040 


— 0,098 


1,7 


0,402 


0,056 


0,169 


0,115 


0,031 


0,086 


1,8 


0,407 


0,041 


0,135 


0,118 


0,023 


0,072 


1.9 


0,411 


0,029 


0,126 i 


0,120 


0,016 


- 0,059 


2,0 


0,413 


0,021 


0,074 


0,121 


0,011 


— 0,041 



Aus diesen Kurven sehen wir bereits die Möglichkeit ein, 
daß Ablösung eintritt. Die Stromlinien enthalten ja die Kurven 
g[^ und Sl^ , welche negativ werden. Also besteht die Wahrschein- 
lichkeit, daß negative Stromlinien auftreten. Diese aber liegen 
gerade inmitten der Grenzschicht. 



T 



\ 
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12. Es gibt noch einen zweiten Weg, auf dem es möglich 
ist die Differentialgleichnngen 

zn integrieren. Differentiiere ich nämlich 2n mal nach iy, so er- 
halte ich 

Die exakte Lösang ist 

Xhirch mehrmalige Integration kann ich dann Tn, Si, L gewinnen. 
Die Integrationskonstanten sind mir bekannt, da ich für iy = 
den Wert jeder Ableitung von f^ bestimmen kann. 




Fig 10.. 

Um nun nicht gezwungen zu sein, bei dieser Eechnung /J, 2n 
mal zu differentiieren, so suche ich das gewonnene Resultat durch 
partielle Integration umzuformen. 

j^«+., ^ e-'''[e'''.4/T-"-/2'?c''*.4/-|.«-'di,], 
/2,/c'''.4/;*»-'>rf,; = [2i,e'''.4/-;,'"-"-/2e'''(l + 2V).4/''"-'di?] 
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/2 (1 + 2V) e^\ Af-^' dn = [-2 (1 + 27/') e''" . 4/';,'»-« -/ (2t; (1 + 2»;') 

g (•»+.) = 4^««-«> _ 2^ . 4/;j— > + 2 (1 + 2ij') . 4/^?-« - 4 (3»; + 2ij») . 4/"»"-« 

-I 1 e J (rationale Funktion von rf) e 4f^ drj. 

Wendet man dann dies für jedes einzelne Glied, in welchem die 
Ableitung von f^ mit einer rationalen Funktion von rj multipliziert 
auftritt, an, so kann man schließlich die gewünschte Ableitung von 
g„ erhalten, ohne /"„ diflFerentiieren zu müssen. Für n = 2 z. B. 
erhielte man: 

g(5) ^ 4/^"-[2i2.4/''-2/V'^^] + [2{l + 2i?').4/^-2/4i?4/'di?] 
= 4r-2i7.4r+4(l + i?').4/--/\(5i? + 2i?»).4/dij 

r = ^r-2nAf+f2Ci + 27i'Wdn^fdnf'-+fff'" u.s.f. 

Man sieht, daß diese Methode besonders für größere Werte 
von n weit komplizierter ist, als das Kuttasche Verfahren. 

13. Wir wenden uns nun der Berechnung von jr, zu, die ganz 
analog der vorhergehenden zu führen ist. %^ genügt der Diffe- 
rentialgleichung : 



dl] 



' ^ dif drj r \ dri [öxötti r dx öiq \ dz öiy* 

ärj [öxÖT] r dx dr] \ dx dif J 

Führen wir für Xo ^^^ Xi ^^® entsprechenden Werte ein: 

du . , ör . 

Xo = rü^o, Xi = ^^0^5'-+^ ö^^>*» 

so geht die Gleichung über in: 



— 30 — 



% + 2n^-8^ = 









+ 



r 



sodaß wir zu setzen haben: 

Dabei hat t^ folgenden Gleichungen zn genügen: 

g;; + 2i?r;. - 8e;, = 4 1 g;. g:- g;' g.. j = 4/-.» 

C+2<c-8gic = 4j-gx.-g.C+3gig;»-g':g..-2g;'g.»} = 4/;. 
c + 2i?g;. - 8g;, = 4 1 gi gl, - r; g.» j . - 4/., 

g;:+2<,-8g;. = 4 {-g.g;»-g.c| = 4/.. 

mit den Graizbedingtmgen : 

für ti = 0: g„ = 0; gj. = 

S.» = 0; g;» = 

5.. = 0; g;. = 
tu = 0; g;. = • 

e.. = 0; g;. = 

für ,, = c»: g;. = 0; g;» = 0; g;, = 0; g;, = 0; t'„ = 0. 

über die numerische Integration ist genaa das gleiche zn 
sagen wie oben. Wir berechnen zunächst ^'^{ri = 0) und J;,(iy = co). 
Diesmal stellen sich die Einflaßkoeffizienten d^ und O*,^ durch fol- 
gende Gleichungen dar: 

K = -A- li^V + 2ij') + (3 + 12^ + 4,?*) e"'/ "«-"'(?,,} 
ö\n ( Ol 
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Die weiteren Rechnungen sollen nicht ausgeführt werden, da sie 
analog den früheren sind. Das Resultat ist folgendes : 

^i(v = 0) = -0,248 g..(i? = oo) = -0,233 



r;i(^ = 0) = -0,068 
C(^ = 0) = -0,029 
f^Liv = 0) = -0,022 
C (^ = 0) = 0,036 



Uil = oo) = -0,043 

tUv = oo) = -0,142 

tuiV = oc) = -0,023 

^^^ (1? - oo) = 0,018. 



Man ersieht schon aus diesen Randwerten, daß der Typus der 
Kurven ganz ähnlich dem sein wird, der schon bei f, auftrat ; nur 
daß bei J^ die Kurven flacher zu werden scheinen. Um die Diffe- 
rentialgleichungen integrieren zu können, brauchen wir noch die 
Zahlenwerte der eingeprägten Kräfte. Dieselben sind in der 
nächsten Tabelle zusammengestellt, in der auch d'^^ und ^,^ zahlen- 
mäßig angegeben sind (Fig. 11): 



f,0 " 



• 1 




Fig. 11. 



n 


4/., 


4/".» 


4/-.. 4/-„ 


4/;. 

1 


*«, 


^.. 



0,1 



0,068 



0,032 



— 0,048 



0,004 



0,011 


1 

0,744 



— 0,026 
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V 


4/.« 


4/. 


4/.C 


^u 


4^. 


K 


*.« 


0,2 


0,244 


0,103 


— 0,154 


0,015 


— 0,046 


- 0,658 


— 0,044 


0,3 


0,494 


0,191 


— 0,236 


0,035 


0,099 


- 0,423 


0,058 


0,4 


0,787 


0,277 


- 0,265 


0,062 


0,166 


— 0,323 


— 0,068 


0,5 


1,093 


0,350 


- 0,206 


0,094 


- 0,235 


— 0,248 


— 0,075 


0,6 


1,386 


0,403 


0,051 


0,128 


— 0,292 


0,193 


0,081 


0,7 


1,645 


0,435 


0,190 


0,160 


— 0,326 


0,151 


-0,085 


0,8 


1,851 


0,445 


0,490 


0,188 


— 0,327 


- 0,119 


-0,088 


0,9 


1,994 


0,437 


0,814 


0,207 


— 0,292 


0,094 


0,091 


1,0 


2,064 


0,412 


1,119 


0,217 


— 0,226 


- 0,075 


— 0,092 


1,1 


2,064 


0.377 


1,384 


0,217 


— 0,137 


0,060 


- 0,094 


1,2 


1,993 


0,335 


1,577 


0,207 


— 0,038 


— 0,049 


-0,095 


1,3 


1,879 


0,288 


1,665 


0,190 


0,060 


0,040 


-0,096 


1,4 


1,719 


0,241 


1,689 


0,167 


0,145 


0,032 


0,097 


1,5 


1,530 


0,198 


1,649 


0,142 


0,210 


-0,027 


-0,097 


1,6 


1,323 


0,156 


1,528 


0,116 


0,251 


— 0,022 


0,098 


1,7 


1,111 


0,121 


1,359 


0,091 


0,269 


0,018 


0,098 


1,8 


0,913 


0,091 


1,181 


0,069 


0,268 


- 0,015 


— 0,099 


1,9 


0,732 


0,067 


0,977 


0,049 


0,253 


0,012 


— 0,099 


2,0 


0,575 


0,049 


0,769 


0,033 


0,230 


0,011 


— 0,099 



«i-- 



.C.2'- 




i,o 



^-»7 



Fig. 12. 



rf-- 




Fig. 13. 
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Sofort können wir nun die Differentialgleichungen integrieren, 
indem wir zuvor 

einmal nacli 1] integrieren: 



Wir wenden wie oben auf die Differentialgleichung das Euttasche 
Verfahren an und berechnen die £ und S'. Das Resultat ist das 
folgende (Fig. 12—16): 




-«»• 



Fig. 14. 




Fig. 15. 



-#,/.- 




Fig. 16. 



fc 



Hf 
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v 


e. 


r 

»Sa 


9ta 


Si* 


bsj 


b!8 


0,1 


0,001 


0,019 


— 0,258 


0,000 


0,007 


- 0,069 


0,2 


0,005 


-0,046 


-0,280 


0,001 


— 0,014 


0,069 


0,3 


- 0,011 


- 0,077 


0,289 


0,003 


0,020 


0,064 


0,4 


0,020 


0,108 


0,271 


0,006 


— 0,026 


0,066 


0,5 


0,032 


0,127 


0,230 


-0,008 


— 0,081 


_ 0,042 


0,6 


0,046 


0,148 


— 0,198 


0,012 


0,035 


0,027 


0,7 


0,061 


— 0,165 


— 0,148 


0,015 


- 0,037 


0,011 


0,8 


— 0,079 


0,177 


0,089 


— 0,019 


0,037 

• 


0,004 


0,9 


— 0,097 


0,188 


0,031 


0,028 


0,036 


0,017 


1,0 


— 0,113 


0,183 


0,281 


- 0,026 


- 0,034 


0,027 


1.1 


0,133 


-0,177 


0,080 


0,029 


0,081 


0,084 


1,2 


0,150 


— 0,167 


0,121 


0,032 


0,027 


0,089 


1,8 


0,166 


0,164 


0,152 


0,035 


— 0,023 


0,089 


1.4 


0,181 


0,137 


0,172 


0,037 


0,019 


0,037 


1,5 


0,194 


0,119 


0,181 


0,039 


— 0,016 


0,032 


1,6 


0,205 


0,101 


0,179 


0,040 


0,018 


0,029 


1,7 


0,214 


0,084 


0,170 


0,041 


0,010 


0,025 


1.8 


— 0,222 


0,068 


0,154 


0,042 


0,008 


0,020 


1,9 


0,228 


0,058 


0,184 


0,043 


0,006 


0,016 


2,0 


0,282 


-0,041 


0,118 


0,043 


0,005 


0,018 



n 


»2« 


t' 

»sc 


»sc 


Ö2d 


t' 


bsd 


0,1 


0,000 


- 0,002 


0,032 


0,000 


0,002 


0,023 


0,2 


0,001 


0,006 


0,044 


0,001 


0,005 


0,026 


0,3 


0,002 


0,012 


0,072 


0,001 


0,008 


0,031 


0,4 


0,003 


0,020 


0,105 


0,002 


0,011 


0,035 


0,5 


0,006 


0,038 


0,186 


0,003 


0,018 


0,036 


0,6 


0,010 


0,047 


- 0,159 


0,005 


— 0,016 


— 0,036 
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V 


92C 


t' 


iie 


ttd 


tu 


t: 


0,7 


0,015 


- 0,065 


— 0,182 


0,007 


— 0,019 


0,032 


0,8 


- 0,023 


— 0,083 


— 0,180 


-0,009 


0,021 


0,026 


0,9 


0,032 


0,100 


0,160 


- 0,011 


0,0-22 


— 0,011 


1,0 


0,043 


0,114 


0,116 


- 0,013 


0,023 


0,008 


1,1 


0,055 


— 0,123 


0,071 


0,015 


0,021 


0,026 


1,2 


0,067 


0,128 


— 0,017 


- 0,018 


— 0,018 


0,037 


1,3 


-0,080 


- 0,126 


0,045 


-0,019 


0,014 


0,041 


1,4 


0,09S 


0,119 


0,094 


0,020 


0,010 


0,036 


1,5 


0,104 


0,108 


0,140 


- 0,021 


-0,007 


0,024 


1,6 


0,114 


0,092 


0,165 


0,021 


— 0,005 


0,017 


1,7 


— 0,122 


0,075 


0,172 


0,022 


0,004 


0,013 


1,8 


0,129 


— 0,068 


0,163 


— 0,022 


0,003 


0,010 


1,9 


— 0,133 


0,044 


0.108 


- 0,022 


0,002 


0,007 


2,0 


0,138 


— 0,037 


0,062 


0,022 


— 0,001 


0,005 



1 

V 


0,1 


0,2 


0,3 0,4 


0,5 


•0,6 


tu 

t',. 

tu 


0,000 

; 0,004 

0,037 


0,001 
0,007 
0,040 


0,002 
0,011 
0,041 


0,003 
0,014 
0,035 


0,004 
0,017 
0,022 


0,006 
0,018 
0,005 



'; 


0,7 0,8 


0,9 


1,0 


1,1 


1,2 


1,3 


tu 

tl» 

tu 


0,008 
0,018 
0,010 


0,010 
0,017 
0,020 


0,011 
0,014 
0,028 


0,013 
0,011 
0,031 


0,013 

0,008 

— 0,028 


0,014 
0,006 
0,019 


0,015 
0,004 
0,014 



n 


1,4 


1,5 


1,6 


1,7 


1,8 


1,9 


2,0 


92« 

t" 
980 


0,016 

0,003 

— 0,011 


0,016 
0,002 
0,007 


0,017 
0,002 
0,005 


0,017 

0,001 

— 0,003 


0,017 
0,001 
0,003 


0,017 
0,000 
0,002 


0,018 
0,000 
0,001 



3* 
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Man sieht, daß die Kurven £, bereits sehr flach werden, und 
es wird daher nicht erforderlich sein, daß wir weitere g aus- 
rechnen. Es soll nur noch der Wert von gj^(iy = 0) bestimmt 
werden; dieser ist wieder leicht durch das Greensche Verfahren 
gefunden. Wir können mit Hülfe dieser Koeffizienten für die Ab- 
lösangszeit and für die Ablösungsstelle eine größere Grenauigkeit 
erreichen. Wir haben zunächst die Differentialgleichungen für f, 
aufzustellen und gehen dazu aus von der G-leichung für x^ Es ist 

äij' ^^ dij' dt) 

r I öi^ \ dxdrj r dx öi; / dx drf 

I h, ( d*x, 1 dr dz, \ dl, d'x. 
dti Xdxdti r dx di] ) dx dt]* 



+ 



h. ( ^\ 1 ^r dx, 
dtj \dxdi] r dx dtj 



\ äxld^xo] 
) dx dt}* J 



Führen wir für Xo> Xu Xt ^^ ^ben gefundenen Werte ein, so ge- 
langen wir durch eine ausgedehntere Rechnung zu einer Differen- 
tialgleichung für Xii die erkennen läßt, daß wir Xt folgendermaßen 
darzustellen haben: 

. .. ö'« ft . -. ö'ö dr , , _, da dV ^ 
,1 „»da (dry , .« ÖV . ^ 1 ^ ar dv . 

+ V"(dx) f- 
Dabei hat g, der Differentialgleichung 

C+2i?g;;-i2ei = 4A 

zu genügen mit den Grenzbedingungen: 

für n =^0: g. == 0; g; = 0, 
für iy = CX5: f^ = 0. 

Zur Berechnung von g, (i^ = 0) ist nur die Kenntnis von 4/*, 
notwendig. Aus der Differentialgleichung für Xt folgen für /*, die 
Ausdrücke : 
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/ I» ©0 ©«• ©0 ©ae • ^^ ©0 ©ac ©O ©!• ^ ©O ©«C ©la ©In ©la ©l6 

/ IC ^== "~ bo b»6 + 2 So £ta + 3 5o Ss6 ■" 2 g^ Sj« "" 2 g^ 5»» ~" SiajCr« + bia ; 
/m ^^^^ '0 ©«i ©oos&i 

/ 8« ^ »0 ©26 • ©0 9tC ©0 ©»6 ©0 ©SC ©la ©iJ • ©la ©ift J 

/ 8/ ^'^^ ©0 bjrf + bo bjc "r 4 bo bfrf "" bo bjc *' ©o bjd H" bia Si6 bi» ©i6 J 

/ 89 ^^ ©0 Öf c ©0 ©2# ©0 ©2« • ^ ©0 ©!• '^ ©0 ©1« ©la © ij ©W ©lÄ 

+ 25IJ — 2g,j5,^; 

/ 8A ^'^ ©0 9td ©0 ©Jd 

tti ^ ©0 bicf — bo birf "^ ^ bo b«« ^ bo b«« bifr b i6 + bi6 

#8* ^^ bo bj« ^ bo hu • ©0 ©8« »i6 • 

Die zelm Funktionen sind mit einem Einflaßkoefflzienten ^,0 za 
multiplizieren und über sie von bis oo zu integrieren, dann haben 
wir die gesuchten g,(?; = 0). Für ^^^ ergibt sich der Ausdruck: 

^30 = ttV j(15 + 90V + 60i^*+8V)e^' P e""^' dri + {3Sfi + 2Sri'+4rj') 



Die folgenden Tabellen enthalten die Werte von 4/", und d' 
(Fig. 17 u. 18). 



so 




*— >7 



«.*.. 



Fig. 17. 
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r.9 




Fig. 18. 



v 


4/-« 


4/^« 


4/-.C 


4:f^ 


4/-.. 


4/-^ 


0,1 


0,040 


0,010 


0,034 


0,002 


0,008 


0,002 


0,2 


0,116 


0,007 


0,092 


— 0,007 


0,025 


0,012 


0,3 


0,183 


0,082 


0.130 

• 


0,014 


0,045 


0,018 


0,4 


0,236 


0,209 


0,146 


— 0,024 


0,061 


0,020 


0,5 


0,24B 


0,387 


0,120 


0,036 


0,063 


0,011 


0,6 


0,212 


0,594 


0,049 


0,047 


0,052 


0,018 


0,7 


0,132 


0,688 


0,046 


0,057 


0,020 


0,070 


0,8 


0,008 


0,725 


— 0,166 


— 0,065 


— 0,028 


0,159 


0,9 


— 0,155 


0,632 


- 0,292 


— 0,069 


0,088 


0,253 


1,0 


0,347 


0,405 


0,410 


— 0,070 


- 0,153 


— 0,357 


1.1 


- 0,528 


0,123 


— 0,505 


-0,069 


— 0,209 


0,450 


1,2 


— 0,726 


0,282 


— 0,574 


— 0,064 


- 0,262 


0,504 
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V 


4/'.. 


4/"« 


4A. 


4/-« 


4/-.. 


4/^., 


1,3 


0,881 


0,725 


0,618 


0,057 


0,291 


- 0,526 


1,4 


1,002 


1,084 


0,603 


0,050 


0,300 


-0,484 


1,6 


1,083 


— 1,434 


0,578 


0,044 


0,284 


- 0,415 


1,6 


1,111 


1,650 


0,540 


0,037 


- 0,262 


— 0,347 


1,7 


1,096 


— 1,746 


-0,482 


0,030 


— 0,225 


— 0,283 


1.8 


— 1,031 


— 1,698 


- 0,411 


0,024 


-0,180 


— 0,222 


1,9 


- 0,939 


- 1,389 


- 0,337 


— 0,019 


-0,140 


0,168 


2,0 


0,823 


— 1,003 


-0,255 


0,015 


-0,102 


0,150 



V 


^U 


4/-.» 


4/^., 


4/^« 


*« 

















1 


0,1 


0,008 


0,000 


0,004 


0,004 


0,700 


0,2 


0,027 


0,002 


0,014 


0,017 


0,495 


0,3 


0,031 


0,005 


0,037 


0,040 


0,353 


0,4 


— 0,003 


0,010 


0,067 


— 0,064 


0,254 


0,5 


0,072 


0,016 


0,102 


— 0,088 


- 0,184 


0,6 


0,193 


0,024 


0,138 


— 0,106 


- 0,135 


0,7 


0,368 


— 0,032 


0,167 


— 0,104 


— 0,099 


0,8 


0,554 


— 0,040 


0,183 


— 0,091 


0,074 


0,9 


0,723 


0,048 


0,179 


-0,064 


— 0,056 


1,0 


0,831 


0,054 


0,149 


— 0,036 


0,042 


1,1 


0,861 


— 0,055 


0,098 


— 0,010 


0,032 


1,2 


0,801 


— 0,046 


0,041 


0,016 


0,024 


1,3 


0,644 


— 0,038 


— 0,000 


0,018 


— 0,018 


1,4 


0,449 


— 0,027 


— 0,016 


0,012 


0,014 


1,5 


0,185 


0,019 


— 0,005 


0,008 


— 0,011 


1,6 


— 0,054 


0,013 


0,000 


0,006 


— 0,009 


1,7 

1 


0,248 


— 0,010 


0,002 


0,004 


0,006 


1,8 


0,308 


0,007 


0,002 


0,002 


— 0,004 


1,9 

1 


- 0,236 


— 0,004 


0,002 


0,001 


— 0,004 


2,0 


0,094 


0,001 


0,000 


0,001 


-0,003 
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Für i'lifj = 0) ergeben sich dann folgende Werte: 

t'^iv = 0) = -0,017 
fUiv = 0) = -0,031 
C(V = 0) = -0,006 
tu (1? = 0) = 0,005 
Civ = 0) = -0,003 
rAv = 0) = 0,008 
r; (1? = 0) = - 0,022 
t'»,iv = 0) = 0,003 
C (ij = 0) = - 0,013 
C (^ = 0) = 0,009. 

14. Über die höheren Glieder der Entwicklung von ^ können 
wir noch kurz folgende Bemerkangen machen; und zwar wollen 
wir znerst untersuchen, wie die Anzahl der 5„ mit wachsendem 
n wächst. Betrachten wir in Xn ^i® Koeffizienten der g„ , so finden 
wir, daß die Dimensionen derselben die folgenden sind: 



für w 


r 


0: 


rtt 


für n 


= 


1: 


X 


für n 


= 


2: 


rü* 


für n 


== 


3: 


rü* 


rü"^^ 








^n 





Also allgemein für n: 

Aus diesen Dimensionen folgt die Anzahl der Glieder daraus, 
daß sie eben so groß sein muß, als man verschiedene Ausdrücke 
von gleicher Dimension büden kann. Dabei ist aber die Ein- 
schränkxmg zu machen, daß x nur als Größe auftreten darf, nach 
der differentiiert wird. Daraus folgt sofort, daß die Anzahl der 
Glieder sich so verteilt: 

für w = 0: 1 Glied: r»; 

für » = 1: 2 Glieder: rtl-r— , ä'-^t-; 

ox dx 

f ür n = 2: 5 Glieder: ^«*(-^), ^«**-tt-, «• s. f. 
für w = 3 : 10 Glieder : ril. (-^j , rü' -^-^ , u. s. f. 
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Diese GliederzaUen ergaben sich schon in nnseren obigen Rech- 
nnngen. Man kann sich nun überlegen, daß X4 ^^^ ^7 Grliedem 
bestehen mnß. Bilden wir die Differenzen der Grliederanzahlen 
für aufeinander folgende n, so ergibt sich das Resultat, daß die 
Anzahl der G-Iieder eine aritbmetische Reihe bilden, deren Dif- 
ferenzen die ungeraden Zahlen sind. Es steigt also die Länge der 
Xn i^t wachsendem n ganz bedentend, und damit ist auch der 
numerischen Berechnung weiterer Grlieder ans praktischen 
Gründen eine Schranke gesetzt. 

Über i'^{rj = 0) und g^(i^= 00) können wir aussagen, daß 
sie mit wachsendem n kleiner werden, daß daher der Einfluß der 
höheren Glieder auf Ablösungsstelle ^) und Ablösungszeit geringer 
wird. Ersteres können wir daraus schließen, daß der Einfluß- 
koeffizient O'^o sich mehr und mehr an die 17 -Achse anschmiegt. 
Andrerseits ist 

I ^^ \ nähert sich also der Asymptote -j — —^ , ohne vorher im 

Endlichen diesen Wert überschreiten zu können. Es muß daher 
die Kurve ^ flacher und flacher werden, damit aber auch ihre 
beiden ersten Ableitungen. 

15. Fassen wir nun kurz unser Resultat bis hierher zu- 
sammen. Wir hatten den Ansatz gemacht 



OD 



^ == 2\ß'Z^fxA^,v) 



woraus : 



y 

2\Joct 



1 OD dy 
r -5^ dri 



folgt. Die Entwicklungen sind bis v =: 2 durchgerechnet, sodaß 
wir folgende Gleichungen haben: 



1) Für die Ablösangs s t e 1 1 e gilt diese Bemerkung nur, falls nicht die in 
Betracht kommende Zeit so groß ist, daß die Potenzen der Zeit, die als Faktoren 
auftreten, überwiegen (vergl. S. 48 u.). 
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+ 1* (rü ( 



"(0S-+''«'^^» + «' 



_. d»» 



, dr du 



dx dx **' 



S, 



ti 



r\drj dt] drj drj ) 



= «e: +< 



+ < 



rlö^j ^- + 7" \äJJ d7^'» + " ^^^^' 



das da;* *•* 



+ « 






., 1 , d*a dr ,, 1 , öS öV 

^^ "•" 7" "ä^ d7 ^»* ■*" 7" öJl«*^ ^" 



-^ ••(£)■ 4 

Die dritte Potenz ist ans folgenden Grründen hinzagefUgt: Die 
Ablösnngsstelle ist gegeben darch 

du 



öy 



= für ij 



Also erhalten wir fnr dieselbe die folgende Gleichnng: 



du 

öy 



= = -^\üC 



2\/xt { 



C (0) + / 



«|^C(0) + ««i£C(0) 









öa;» 



öwc" öj; 
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, 1 ,Ö0 /drV 

+ — «ä- 
r ax 






Dividieren wir durch 



ü 



2\ß 
neten Werte ein, so erhalten wir: 

= 1,128 + /[l,614|5 + 0,169« ^ ^'" 



und führen wir für gj^(0) die berech- 



dx 



r dx 



- f [0,248 (g)' + 0,068 ü S + 0,029 ö 1 ^« ^»^ 



r ö« d» 



öa;' 
+ 0,022 .•i;^-O,OB6J,»-(0] 

-. [0,017 ©V 0,031 «i(£)'i^ + 0,«»a 



>x ^. ö*« 



, 1 d*a dr 



da d'ä 
dx dx* 

., 1 da ä'r 



-0,005«'^ + 0.003«' -^^-0,008«'^^^ 



+ 0,022 4 «' 1^ (P! - 0,003 i «• -S- + 0,013 4 « 



-0.010^ «•(£)] -[.•]. 



Dies ist die endgültige Ablösangsgleichnng für be- 
liebige Rotationskörper. 

16. Ist mir die Stelle bekannt, an der die Ablösung znerst 
eintritt, so kann ich mir aas dieser Gleichung die „Ablösangszeit" 

berechnen. Dies würde gelten für Körper, die weder in Spitze 
noch Kegel auslaufen und deren Ende konvex ist ; denn von diesen 
Körpern können wir aussagen, daß die Ablösung zuerst am hin- 
tersten Ende statthat, falls nicht an einer weiter nach vorn lie- 
sehr groß ist. Bei den hier ausgeschlossenen 



genden Stelle 



dx 



Körpern kann es sehr wohl geschehen, daß die Ablösung an irgend 
einer anderen Stelle eintritt. In diesem Falle kann man die Ab- 
lösungszeit so bestimmen, daß man für verschiedene Ablösungs- 
stellen die Zeiten' berechnet. Trage ich mir dann die Zeit als 
Funktion dieser auf, so muß die Kurve an der Stelle, wo die Ab- 
lösung zuerst eintritt ein Minimum haben. 
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III. Anwendung der Resultate auf die Kugel. 



1. Für die Eagel ist 



Ä = J Fsin 



X 



R' 



wenn V die Strömungsgeschwindigkeit der Flüsssigkeit im un- 
gestörten Zustande oder auch die Geschwindigkeit, mit der die 
Kugel durch die ruhende Flüssigkeit bewegt wird, ist. Mit R 

bezeichnen wir den Kugel- 
radius, mit X die ßogen- 




Fig. 19. 



also: 



da 
dx 

dx" 
ö^ 



3^J^ x_ dr 

2 R^""^ r'^ ~d^ 

2 R' ^^^ R ' dx' 

IZ ^ ^ 

2 R' ^^^ R ' dx' 



länge längs eines 
Kreises (Fig. 19). 


größten 
Dann 


ist 








r 


== Jisin-ö 


• 

'j 




X 

cos^ 






1 
R 


. X 

srn^ 






1 
R* 


cos^ 


• 





2. Unsere Ablösungsgleichung lautet, wenn wir durch 1,128 
dividieren : 

xlS r 



1 + 1,423 cos ^ + 0,150 cos |-] J ^ ^ 



X 



0,220 cos' -^ 



0,060 sin' ^ + 0,026 cos« -^ - 0,020 sin' 4 



- 0,032 cos' 



r\\2R') 



0,016 cos» ^ - 0,028 cos' ^ - 0,006 sin* -^ cos ~ 

li li Jti li 



+ 0,005 sin' I- cos |r - 0,002 sin' l cos ^ + 0,007 sin* ^ cos -^ 

K Xt li U RR 
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+ 0,020 cos* ^ + 0,002 sin' ^ cos -^ - 0,011 sin' -^ cos ^ 

M UM XV li 



-0,008co8-|](|^/y = 0. 



Der Einfachheit halber führen wir neue Variable ein, indem wir 
setzen : 

X _ ^ 3 V _ 
R ~ ^' 2R^ - ^' 

Bei der Kngel tritt die erste Ablösang sicher bei 

X = n 

ein, sodaß wir für die Ablösnngszeit [TJ die Grleichang erhalten: 

1- 1,573 [2] -0,214 [27 ~0,(X)1[T]' == 0. 

Für [T\ ergibt sich daraus 

[T] = 0,589 
also 

[t] = 0,392 ^ ■ 

Man sieht, daß die dritte Potenz von [T] für die Ablösungszeit 
nur wenig in Betracht kommt. 

Die Kugel hat, wenn die Ablösung eintritt, bereits einen 
Weg zurückgelegt, der annähernd f des Kngelradius lang ist ; dies 
gilt für jede Greschwindigkeit und jeden Radius. 

3. Wir berechnen uns nun noch die Stromfunktion ^. Es 
war 

Auf die Kugel angewendet: 
n, = 2V^[||^sm'|-f. + | >^ sin' |- cos |- < (6.. + e J 



+ 



^ -^ sin' |- 1* |co8'|- (g,. + 5,. + o - sin' |-;(g;» + u}] 



Setzen wir wieder 
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— — Y- — — t — T 



so veremfacht sich die Gleiclmng dahin: 



+ r { cos* X (S„ + e.. + gj - sin« X (U + f«,) { ]. 
Endlich führen wir auch für ^ eine neae Variable ein: 



sodaß wir schließlich die einfache Gleichung haben: 

y = V^sin* X K. + r cos X {t.a + 5J + ^ I cos» X (g,„ + g.. + S,.) 

-8in'X(g.» + g.,)|- 

Wir können nun leicht für T irgend welche Werte einsetzen und 
damit ?F als Funktion von tj berechnen. Aus dieser Kurvenschaar 
findet man auf einfache Art die Stromlinien, d. h. ^ als Funktion x. 
4. Es empfiehlt sich, um genauer den Vorgang der Ablösung 
zu studieren, zunächst zuzusehen, wie die Dicke der Wirbelschicht 
mit der Zeit wächst. Wir betrachten zu dem Ende die Strom- 
linie 5** = 0, und fragen uns hier wiederum zuerst, wie bewegt 
sich der Ablösungskreis auf der Kugel fort. Wir setzen in der 
Ablösungsgleichung für X die Werte von 7t an abwärts ein und 
erhalten so für T resp. t Gleichungen dritten Grades, die man 
mit Hülfe des Newtonschen Interpolationsverfahrens unschwer 
lösen kann. Für die Ablösungsgleichung schreiben wir: 

1 + 1,573 cos XT- (0,214 cos' X - 0,080 sin» X) T* 

+ (0,001 cos* X - 0,005 sin* X cos X) T« = 

oder in einer zur Rechnung geeigneteren Form : 

1 + 1,573 cos XT~ (0,294 cos« X - 0,080) T' 

- (0,003 cos' X- 0,005) cos XT» = 0. 

Führen wir die Rechnung für verschiedene Werte von X aus, 
so erhalten wir folgende Kurve (Fig. 20): 



X 

1 


ISO» 


175» 


170« 


165» 


160» 


155» 


160» 


T 

4 


0,589 
0,392 


0,591 
0,394 


0,598 
0,399 


0,610 
0,407 


0,628 
0,419 


0,653 
0,435 


0,685 
0,457 
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X 


14B» 


140« 


135» 120« 


IIB»' 


110« 


90« 


T 


0,727 


0,782 


0,854 


1,280 


1,599 


2,258 


OO 


'l 


0,484 


0,521 


0,569 


0,853 


1,066 


1,505 


oo 


* \ \ 








\ 


\ 
\ 
\ 














.. t»' 



* 4» 



9.S9 



I 



»y 



Fig. 20. 



Die Werte sagen über das Portschreiten der Ablösung fol- 
gendes ans: nach Verlauf der kritischen Zeit 

beginnt die Ablösung plötzlich, und die Ablösungsstelle bewegt 
sich anfangs mit unendlicher großer Geschwindigkeit fort, läuft 
dann langsamer und langsamer und erreicht nach unendlicher langer 

Zeit die Stelle X = -^. Über diesen Punkt hinaus kann sie 

nicht fortschreiten, solange die Dicke der Wirbelschicht hinter der 
Kugel — so wie hier angenommen — noch klein gegen den Kugel- 
radius ist. Die Ablösung tritt zuerst ein, wenn die Kugel sich 
ungefähr um f ihres Radius fortbewegt hat und ist bereits um 
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38^ weitergerückt, wenn die Kugel tun die halbe Badiuslänge, 
also nur nm ein V^o weiter als vorher, sich verschoben hat. Dann 
geht es langsamer , sodaß auf eine weitere Weglänge von einem 
halben Kngelradias nur noch 28^ kommen. In Figar 21 sind die 
vier genannten Lagen der Kugel mit ihren entsprechenden Ab- 
lösungsstellen schematisch wiedergegeben mit Innehaltung der Ver- 
hältnisse (die Dicke der Wirbelschicht ist überhöht). 

5. Die weitere Untersuchung bezieht sich darauf, wie dick 
zu bestinunten Zeiten an den verschiedenen Stellen der Kugel die 
Wirbelschicht ist, d. h. wir untersuchen für verschiedene X den 
Abstand y der Stromlinie !F = als Funktion von t Zur Be- 
rechnung wird die Gleichung der Stromlinie: 

W = sin« X VT [g, + T cos X (S,« + W 

+ T \ cos« X (e.. + f., + gj - sin" X (g,, + U! 

verwendet. Für y = können wir schreiben : 

= to+^cosX(g.. + gJ + r* {cos'X(t,. + 5„ + g„ + C., + U 




Fig. 21. 

Dabei ist g als Funktion von ?; berechnet, und es ist deshalb zu 
berücksichtigen, daß: 

2^2 2 kR 

^ - 2 R^ 

ist. Für X setzen wir die verschiedenen Werte von n bis -5- ein 
und erhalten so eine Kurvenschaar 
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y = /"(X, 0, 

wo X die Rolle eines Parameters spielt (Fig. 22). Die Kurven, 
deren Punkte zahlenmäßig anzugeben wenig Wert hat, sind an- 



nähernd gradlinig, d. h. 



dy 
dt 



rsj 



const., also wächst die Dicke mit 



konstanter Geschwindigkeit. 

Dieses Eesnltat stimmt o£Penbar mit den Erscheinungen nicht 
überein. Man beobachtet nämlich, daß das Anwachsen der Wirbel- 
schicht immer rascher nnd rascher vor sich geht, so daß ihre Dicke 
in kurzer Zeit mit dem Kagelradias vergleichbar wird. Das 
Gradlinige kommt nur daher, daß die Reihe hinter dem Glied T* 
abgebrochen wird; und es ist vorauszusehen, daß für größere T 
die höheren Glieder nicht gegenüber den niederen vernachlässigt 
werden dürfen. Die weiteren Aassagen sind daher nur unter der 
Voraussetzung, daß die Glieder von höherer als zweiter Potenz in 
T vernachlässigt werden dürfen, zu verstehen. 



■Ä 



.»" 



».f 









«i 




Fig. 22. (Die Kurven sind von 5° zu 5® angegeben von 110° — 180**). 

Die Geschwindigkeit ist für X = ä am größten und nimmt 

ic 

mit X ab, bis sie bei X = -^ den Wert Null erreicht. Um einen 

Einblick in die relative Dicke der Grenzschichten zu erhalten, 

Ji 

wollen wir einen speziellen Fall betrachten, t sei gleich 0,5 — , 

dann ist die Ablösung ungefähr um 38^ fortgeschritten, und y ist 

4 



N 
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für X = Ä un 



cm 



gefähr 0,59 \/- 



V 



Nehmen wir 22 = 10 cm und 



y = 10 , dann wird 



sec 



und, wenn wir 



t = 0,5 sec 



X = 0,01 



cm 
sec 



(Wasser) 



annehmen, 

y = 0,059 cm. 

Also verhält sich der Radias zur Dicke des Wirbels wie 1000 : 6. 
Die Störung durch die Reibung erstreckt sich auf ungefähr 0,3 cm, 
und die Verdickung (s. u. 6.) beträgt ungefähr 0,19 cm. Man sieht 
aus diesen Angaben, daß wir die Grrenzschicht bildlich nur dar- 
stellen können, indem wir für y einen Maßstab wählen, der groß ist 
gegen den der X-Koordinate. 

Aus der Kurvenschaar Fig. 22 können wir sofort eine neue 
finden, die uns die Lage der Stromlinie ^ = für verschiedene 
Zeiten angibt. Aus derselben können wir uns ein Büd über den 
Verlauf der Stromlinien ^ ^ machen ; wir sehen zugleich, wie 
sich die Wirbelschicht weiter und weiter hinter der Kugel aus- 
dehnt (Fig. 23). 




f 
Fig. 23. (Die Kunren sind von 0,1 zu 0,1 angegeben). 

6. Die Übersicht über die Wirkung der Grenzschicht wird 
noch weiterhin sehr erhöht, wenn man noch die scheinbare Ver- 
dickung des Körpers zeichnet. Da außerhalb der Grenzschicht 
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Potentialströmimg herrschen soll, so muß es möglich seinj, den 
Körper za finden, der zu dieser reinen Potentialströmung gehört. 
Die in dieser Arbeit berechneten Geschwindigkeitsprofile (u als 
Funktion von y bei konstantem X) näheren sich einer Asymptote 
ft. Diese Asymptote stellt die Greschwindigkeit bei vollständiger 
Potentialströmung dar (Fig. 24). Deshalb muß die Grenze des 
Hilfskörpers an der Stelle y = y liegen, welche (wegen der Kon- 
tinuität) gleichen Inhalt der in Fig. 24 schraffierten Flächen liefert. 
Dieses y entspricht dem Schnittpunkt der Asymptote der Kurve 
tlf als Funktion von y für konstantes X. Die so gefundenen Ver- 
dickungen sind in Fig. 23 eingetragen. 

7. Der eingeschlagene Weg der Untersuchung 1) der Aus- 
dehnung, 2) der Dicke des Wirbels und 3) der Verdickung ist dem 
vorgezogen, auf dem man von den Strömungs Widern selbst aus- 
geht, obgleich der letztere Weg wohl schneller zum Ziele zu führen 
scheint. Bei letzteren Betrachtungen würden die Geschwindigkeiten, 
mit denen sich die Grenzschicht ändert, nicht so klar hervortreten, 
wie sie es mir bei dem oben genannten Verfahren zu tun scheinen. 




Fig. 24. 

8. Dennoch ist es von großem Nutzen, wenn wir für eine 
bestimmte Zeit, in der der Wirbel schon weiter ausgebildet ist, das 
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Bild der Stromlinien konstmieren und daran die Greschwindigkeiten 
studieren. In Fig. 25 ist dies geschehen für die Zeit t = 0,6 sec 
bei einer Geschwindigkeit von der Kadinslänge pro Sekunde. Die 
Überhöhung der Maße in der y- Richtung ist ungefähr 30 fach. 
Für die Stromlinien erhalten wir die Gleichimg: W = const., wo 

W = 0,775 sin« X [g, + T0,6 cos X (g,. + U 

+0,36 {cos* x(g„+g,,+g,,+s,,+ej-(e,,+gj|]. 

Damit haben wir ^ als Funktion von rj und damit von y bei 
konstantem X, Das Kriterium dafür, daß Wirbelbildung eintritt 




Verdickung 



y = 



X = « 



^-f 



Fig. 25. 



1 
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ist das, daß V{rjD negativ wird. Trägt man aaf einer Anzahl Radien 
die Fonktionswerte der VP aaf, so läßt sich hieraus sofort y als 
Funktion von X für jede Stromlinie ^ «= const. entnehmen. 
"Wenn wir noch für die verschiedenen X die Geschwindigkeits- 
profile auftragen (diese sind in Fig. 26 zur Darstellung gebracht), 
so sieht man, daß in dem Wirbel vorerst nur relativ kleine Gre- 
schwindigkeiten auftreten. Außerhalb desselben steigt die Gre- 
schwindigkeit sehr schnell auf den Fotentialwert an. Man kann 
auch sofort die Wirbelstärke an den einzelnen Punkten aus der 
Figur entnehmen; denn es ist 

, dv du 

rot ü = 



dx dy 
Nun ist -^ — r\f c, kann also vernachlässigt werden, sodaß 

rot ö = — 



öy 



Der Wirbel ist gleich dem negativen Gradienten der Geschwindig- 
keit in Richtung der Flächennormale. Die Rotation ist, wie 
man sieht, absolut am stärksten außerhalb der Stromlinie ^ = 
und ist in der Nähe der Kugeloberfläche nur sehr gering. Schließ- 
lich ist in Fig. 25 auch die Verdickung mit eingezeichnet. 



